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Abstract 

o _ 

CN \ We analyze neutral Higgs boson decays into squark pairs in the minimal super- 

CD ' symmetric extension of the Standard Model and improve previous analyses. In par- 

ticular the treatment of potentially large higher-order corrections originating from 



the soft SUSY breaking parameters A\> % the trilinear Higgs coupling to sbottoms, 
and n, the Higgsino mass parameter, is investigated. The remaining theoretical 
uncertainties including the SUSY-QCD corrections are analyzed quantitatively. 
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' ■ 1 Introduction 

P+ 

1) ■ 

^ | The search for Higgs bosons plays one of the most important roles at high-energy collider 

experiments at the Tevatron and the LHC. The Higgs boson is the remnant of electroweak 
symmetry breaking in the scalar Higgs sector of the Standard Model (SM) and its su- 
^ , persymmetric extensions. The minimal supersymmetric extension of the SM (MSSM) 

requires the introduction of two Higgs doublets in order to preserve supersymmetry pp. 
This leads to the existence of five elementary Higgs particles, two CP-even (h,H), one 
CP-odd (A) and two charged (H^ states. At lowest order (LO) all couplings and masses 
of the MSSM Higgs sector are described by two independent input parameters, which are 
usually chosen as tg/3 = v 2 /vi, the ratio of the two vacuum expectation values u li2 , and 
the pseudoscalar Higgs-boson mass At LO, the light scalar Higgs mass has to 
be smaller than the Z-boson mass Mz- Including the one- loop and dominant two- loop 
corrections the upper bound is increased to M h < 135 GeV [2]. Recent first three- loop 
results confirm this upper bound within less than 1 GeV [3] . The couplings of the various 
neutral Higgs bosons to fermions and gauge bosons depend on the mixing angles a and 
p. Normalized to the SM Higgs couplings, they are listed in Table [TJ For large values of 
tg/3 the down-type Yukawa couplings are strongly enhanced, while the up-type Yukawa 
couplings are suppressed, unless the light (heavy) scalar Higgs mass ranges at its upper 
(lower) bound, where the couplings become SM-like. 

The negative direct searches for the MSSM Higgs bosons at LEP2 yield lower bounds 
of M KH > 92.8 GeV and M A > 93.4 GeV. For a SUSY scale M SUSY = 1 TeV the range 
0.7 < tg/3 < 2.0 in the MSSM is excluded by the Higgs searches at the LEP2 experiments 
[I] . Presently and in the future, Higgs bosons can be searched for at the Fermilab Tevatron 
[5], a proton-antiproton collider with a center-of-mass energy of 1.96 TeV, and the proton- 
proton Large Hadron Collider (LHC) with up to 14 TeV center-of-mass energy [5J. 
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cos a/ sin [3 — sin a/ cos (3 sin(/3 — a) 
sina/sin/3 cosa/cos/3 cos(/3 — a) 
l/tg/3 tg/3 



Table 1: Higgs couplings in the MSSM to fermions and gauge bosons [V = W, Z] relative 
to SM couplings. The subscripts u, d denote up- and down-type fermions. 



The MSSM Higgs bosons couple to squarks, too. If the Higgs masses are large enough 
to allow for decays into third-generation squark-antisquark pairs, these decay modes ac- 
quire sizable branching ratios in many MSSM scenarios [7]. Thus, for a reliable theoretical 
treatment of these supersymmetric Higgs decay modes, higher order corrections have to 
be computed and included appropriately In the past the full SUSY-QCD corrections [8] 
and the full SUSY-electroweak corrections [9] to the Higgs decays into squarks have been 
calculated. In particular regions of the MSSM parameter space the corrections turned 
out to be so large that a reliable prediction was not possible without further refinements. 
A first attempt to solve this problem has been undertaken in Ref. [10] by starting from 
a more consistent treatment of the squark masses and couplings at next-to-leading order 
(NLO). However, the approach of Ref. [TD] does not provide a treatment of the squark 
mixing starting from the running soft SUSY-breaking parameters of the squark sector 
only. The topic of this work is a complete and consistent determination of the MSSM 
squark sector and the Higgs couplings to squarks from running MS input parameters in 
the context of SUSY-QCD corrections^. The input parameters can be obtained from the 
renormalization group equations in the framework of a certain SUSY-breaking mechanism. 
Our treatment can also be used for a consistent global fit to supersymmetric observables 
which include observables based on the stop and sbottom sectors of the MSSM. 

This paper is organized as follows. Section [2] describes the systematic determination of 
the squark sector at NLO, while in Section [3] we summarize the SUSY-QCD corrections 
to the MSSM Higgs decays into squarks. Numerical results are presented in Section HI 
Finally in Section [5] we conclude. 

2 Squark masses and couplings 

In this section, we describe in detail the determination of the stop and sbottom masses and 
their couplings to the MSSM Higgs bosons at LO, and extend the setup to NLO consis- 

1 We have chosen the MS scheme for convenience while our analysis could easily be translated to the 
DR scheme as used in spectrum generators. The MS parameters can be obtained from the D R ones by 
simple relations [llj and vice versa. 
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tently. The NLO expressions will be derived for soft supersymmetry breaking parameters 
given in the MS scheme. 



2.1 Sfermion masses and couplings at LO 

The scalar partners Jl,r of the left- and right-handed fermion components mix with each 
other. The mass eigenstates of the sfermions / are related to the current eigenstates 
Il,r by mixing angles 9 f , 



fi = f L cos9 f + f R sm9 f 
h = -/l sin 6>/ + f R cos9 f , 



(1) 



M, 



(2) 



which are proportional to the masses of the ordinary fermions, see Eq.(|3J). Thus mixing 
effects are only important for the third-generation sfermions t, b, f , the mass matrix of 
which is given by 

M 2 h +m) m f (A f - (jtr f ) 
m f (Af-fir f ) Mj R +mj 

with the parameters r& = r T = l/r t = tg/3. The parameter Af denotes the trilinear 
sfermion coupling of the soft supersymmetry breaking part of the Lagrangian, while \x is 
the Higgsino mass parameter and rrif the corresponding fermion mass. The D-terms have 
been absorbed in the parameters r , 



Mj 

JL/R 

Df 

JL 

^Ir 



Mj +D f 

JL/R JL/R 



where M 



Il/r 



M\(l{ h - e f sin 2 9 W ) cos 2/3 

Mfe/sin 2 6Vcos2/3, (3) 
denotes the sfermion masses of the soft supersymmetry breaking part of the 



Lagrangian. Consequently the mixing angles acquire the form 

2m f (A f - fir f ) 



sin 29 



f 



2 
h 



2 
h 



cos 29 



f 



Mj -Mj 

JL JR 
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h 
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J2 



and the masses of the squark eigenstates are given by 

1 



/l,2 



m 2 f + 



JR. 



Mf L + M% T \/{Ml - M}J + Am){A f - pr,) 



Jl Jr 

In the current eigenstate basis the neutral Higgs couplings to sfermions read a^| 

m)gf + M 2 z (hf-e f sin 2 9 w )gt 
m 2 f gf + M 2 z e f sin 2 w g£ 

A f gf) 



(4) 



(5) 



SfLfL 
9 f R fR 
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SlSr 



(6) 



2 In our notation the first index of the neutral Higgs couplings to sfermions defines the incoming and 
the second index the outgoing sfermion at the corresponding vertex. 



with the couplings gf 



1, . . . , 4) listed in Table [2j For the scalar Higgs bosons h, H 



the couplings to sfermions are symmetric, i.e. g 



h H h H 

7 7 — Qi i ■ For the pseudoscalar Higgs 
boson A the diagonal couplings gf - and gf ~ vanish, while the off-diagonal couplings 

JLjL JRjR 

are antisymmetric, i.e. g A - 1 



JrJl 

eigenstates /i^ are given by 



g'~ e . The corresponding couplings to the sfermion mass 

J L J Ft 



h,H 
9f f 

h,H 

^hh 

h,H 

97 7 

fih 

9ff 
^ fill 

^ fifi 



n.n 2 n i h,H ■ 2 n n.a ■ <-»/) 

9i k cos 9 f + 9f Rf ~ R sm f + g f ' T K sm 20, 

h,H 



h,H 

J f L fL 

h,H „■ 2 



h,H 

fhfR 

h,H 



n-.ti -In i n,a 2 n not ■ nn 

9 kk Sm d f + 9 fRfR B * ~ g kfR Sm 29 f 



h,H / ft,.H ft,_H \ • o/l i ftJJ o/l 

= 2 {9 h.fR - g hh ) Sm2 ^ + g hfR COS2 ^ 



g hh g hf~R 



(7) 



For a consistent NLO calculation the relations for the sfermion masses, mixing angles and 
couplings have to be extended to NLO consistently. In the following we will concentrate 
on the stop and sbottom sectors at NLO SUSY-QCD. 
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Table 2: Coefficients of the neutral MSSM Higgs couplings to sfermion pairs. The symbols 
u, d denote up- and down-type sfermions. 



2.2 Stops and Sbottoms at NLO 

The soft supersymmetry-breaking parameters Mg LR (Q ) and A q (Q ) will be introduced 
as MS parameters at an input scale Qq which will in general be of the order of the SUSY 
scale. 

In order to resum large corrections for large values of tg/3 we define the sbottom masses 
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and couplings to the Higgs bosons in terms of the effective bottom masa^] [10J 121 113] 



rh b (Q) 



m b (Q) 
1 + A, 



b 



Ab = 9L ^M §/ /tg/3 /(ra^raf^Mf) 

a , , c 
ab log — + be log - + ca log - 

(a — b)(b — c){a — c) 



with Cf = 4/3, where m b (Q) denotes the running bottom mass in the MS scheme at 
the scale Q and rag , mg 2 , M g are the sbottom and gluino pole masses respectively. The 
effective top mass is identified with the running MS top mass 

m t {Q) = m{Q). (9) 

The use of these effective masses corresponds to the consistent inclusion of the resummed 
and RG-improved Yukawa couplings of the Higgs decays 0° — > ti, bb everywhere in the 
squark sectors, tocQ. The stop/sbottom mass matrix at LO is then given by (q = t, b) 



m\ l {Q q ) + m 2 (Qo) rh q (Q )[A q (Q ) - /ir q ] 



^2 



rh q (Qo)[A q (Q ) - fxr q ] M-(Q ) + m 2 q (Q 



(10) 



where A q (Q ) denotes the running trilinear MS coupling at the scale Q and r b = l/r t = 
tg/3. Analogous to Eqs. ([3]) the D-terms D qL/R have been absorbed in the parameters 

^ L/R (Qo) = M 2 qL/R (Qo) + D qL/R . (11) 

The stop/sbottom mass matrix is modified by higher-order corrections in the diagonal 
and off-diagonal entries. While the corrections to the off-diagonal entries will be treated 
by the renormalization of the mixing angles, we compensate the radiative corrections to 
the diagonal matrix elements by shifts in the soft mass parameters M qL/R (Q Q ), 



~ 2 9 



Mi , (Qo) = M~„ (Qo) + AM, , Mi (Q ) = M~ n (Q ) + AM, (12) 



3 Note that this definition of the effective bottom quark mass differs from the DR definition which 
has been used in Ref. [10] . Our effective bottom mass runs with five active flavours consistently, i.e. all 
heavier particles are decoupled, while the DR definition requires running with the contributions of all 
supersymmetric particles. The scale of the strong coupling a s in within the effective bottom mass 
is identified with the average mass of the sbottom and gluino states in order to account for the NNLO 
corrections to a large extent [14] and is thus different from the scale Q. 

A similar approach has meanwhile been pursued in Ref. [15j for the Higgsino couplings to fermions 
and sfcrmions, too. 
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in order to arrive at simple expressions at NLO for the stop/sbottom masses. Using the 
squark pole masses Tn qi . 2 the tree-level definition 9 q of the mixing angle is given by 



sin 29 „ 



2m q (Q )[A q (Q ) - fj,r q ] 

i i 
m~ — m% 

91 92 



cos 29 n 



Ml(Q )-Ml(Q, 
ml - ml 



(13) 



where the shifted squark mass parameters M qL/R (Q ) have been used with the radiatively 
corrected squark pole masses. This definition will be used as the tree-level-like mixing 
angle with loop-corrected squark masses m qi/2 at NLO, too, but will only play the role of 
an auxiliary parameter in our analysis as will be explained later. 
At NLO the masses of the stop/sbottom eigenstates are given by 



m 



Ql/2 



Am 2 , 

Ql/2 



rh q (Qo) + 2 



mUQo) + M 9 (Qo) 



~ 2 



~ 2 



T\ [M^(Q ) - M qR (Q )] 2 + Am 2 (Q )[A q (Q ) - f,r q ] 



+ Am l /2 



'n/22(^i ) + 5m 2 ; 

ii/za\ q x / 2 J 9l/2 



(14) 



The diagonal parts £11/22 of the stop/sbottom self-energies can be calculated from the 
diagrams in Fig. (TJ 

a 

9 9 



Qi 



Si + 



Qi 



+ 



Qi 



q q 

Figure 1: One-loop contributions to the squark self- energies. 



q.i 



£1 



1/22 ^- i/2 



Cf ~Z + C ° s2 2 ^)^oK- 1/2 ) - sin 2 29 q A {m qyi ] 
+2A (M- g ) + 2A [m q (Q )]+Am 2 i/2 B (m 2 i/2 ;0,m qi/2 ) 



+ 2 



(15) 



Ml + m 2 q (Q ) - m 2 qi/2 T 2M g m q (Q ) sin 20, 5 K~ 1/2 ; M g ,m q (Q )} 

The one-loop integrals are defined as (16] 

d n k -i(47r) 2 /i 2e 



Anlm) 



B (p 2 ;mi,m 2 ) 



(2vr) n k 2 - m 2 
d n k -i(4n) 2 fi 2e 
(2vr)™ [k 2 - m 2 }[(k + p) 2 - m 2 } 



B 1 (p 2 ;m 1 ,m 2 ) = {^o(™i) - A (m 2 ) - {p 2 + m\ - ml)B (p 2 ; m h m 2 )} . (16) 
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The parameter fi denotes the 't Hooft mass of dimensional regularization. The counter 

91,2 



terms bm\ are given by 



5m| /2 = 2m q (Q )8m q + - {sm\ l + 5M\ R ± [{5M\ l - 5M\ R ) cos 20, 



I t%t + TTT^ ) « - ml) sin 2 2^ 



= -C F ^T(1 + e)(47r) e |i + log^j {M| T M 9 m 9 (Q ) sin 20,)} 

+^^y {2m 2 (Qo) =F ~« - <) sin 2 2^} , (17) 

where the tree-level-like mixing angle 9 q of Eq. (1131) has been used and the parameters 
M-: and A a have been renormalized in the MS scheme, 

r.2 



SA q = C F ^T{1 + 6){A<KyM g \- + \og^\ . (18) 



7T ^ ' ^ U ~°Q 2 

The counter term of the effective quark mass m q (Q) for Q = Q is given byEI 
() "'" - -C^r(l + e)(47r)4|-+logC + wl + A, 



^,,,7 ^ [m 2 (Q); Mg, m^] - £ [m 2 (Q); M ff , 
+2M 9 (A g - M r 9 ) _ \ , (19) 



where Ssusy = 1/3 is a finite counter term required to restore the supersymmetric re- 
lations between the Higgs boson couplings to top/bottom quarks and stops/sbottoms 
within dimensional regularization [TTJ. The term A q denotes the correction A& of Eq. (JS} 
for the bottom mass, while it vanishes for the top mass case, i.e. A t = 0. The term 
A;, in the counter term of the bottom mass cancels the leading term for large values 
of tg/3 in the last line of Eq. f fl9|) so that this counter term is free of large corrections. 
Using the NLO corrected squark pole masses of Eq. ( !T4|) and the tree-level-like mixing 
angle 9 q of Eq. ( fl3l) . the shifted squared soft SUSY-breaking squark mass parameters 



3 The quark mass counter term is different from the DR expression used in Ref. 10 . 
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~ 2 



M| (Qo) = My L . (Qo) + AM? can be obtained from the relation 



&1 



M\ L (Qo) = M| £ (Qo) + D 9l = ml cos 2 9 q + m? 2 sin 2 ^ - m\ (Q ) 

M 2 fl (Q ) = M 2 R (Q ) + ^ = m 2 i sin 2 ^ + m 2 2 cos 2 ^-mJ(Qo). (20) 

The tree-level definition of the mixing angle 9 q in Eq. (|T3|) corresponds to the following 
counter term at NLO, 



tg 29 q f <5m g (L4 g gm| - gm| 

' 1 m 9 (Q ) + A 9 (Qo)-/ir g m| - m| ^ 



<K /2 = - s n/ 22 K /2 )- (22) 

However, in order to avoid artificial singularities in physical observables for stop/sbottom 
masses m qi 2 close to each other, in our calculation the mixing angle of the squark fields 
has been renormalized via the anti-Hermitian counter ternfli^, 

l 3fteSi2K) + ^E 12 (m|,) 

OVq = 2 2 ' y 26 ) 

where S 12 denotes the off-diagonal part of the stop/sbottom self-energy (see Fig. [Q de- 
scribing transitions from the first to the second mass eigenstate or vice versa, 

S 12 (m 2 ) = -C F —{M- g m q (Q )B [m 2 -M- a ,m q (Q )] 

7T 

sin 29 I 
+^T A [ A °( m 9~ 2 ) - ^oK)l f cos 2# 5 . (24) 

This implies a finite shift A9 q to the mixing angle 9 q of Eq. ( IT3|) . 

9 q = 9 q + A6 q , A6 q = 56 q -de q (25) 

which modifies the relations of Eq. fl2"U|) by replacing 9 q = 9 q — A9 q . The scale and scheme 
dependence of the input parameters determining the squark pole masses is explicitly 



6 These equations differ from the corresponding inconsistent relations in terms of the on-shell definition 
of the mixing angle used in [lOj . Our approach starts from the consistent tree-level relations and incor- 
porates the anti-Hermitian definition of the mixing angle in a finite shift from the tree-level value 8 q . An 
alternative and equivalent option would be to compensate the shift A9 q by a shifted trilinear coupling A q . 
Moreover, it should be noted that the shifted parameters M^ l (Qq) and Mj^Qo) are not equal to each 
other any more in accordance with the different radiative corrections to the stop and sbottom masses. 

7 Alternatively the mixing angle can be renormalized by the counter term S9 q — Sfte£i2(Q 2 )/(TO? 2 — m^) 
with Q 2 = or Q 2 = m 2 ~ 2 [T7| which leads to a result free of the artificial singularities for stop/sbottom 
masses m qi 2 close to each other, too. Another option free of these singularities is provided by taking the 
residue of Eq. ([25)) for equal masses for the mixing angle counter term [T5] . 



8 



compensated by the corrections Am| 12 to the tree- level relations of Eq. (JHJ), while the 

scale dependence of the mixing angle is compensated by the finite shift A6 q in Eq. ( )25l) . 
It should be noted that the corrected mixing angle fulfills the relation sin 2 26 q + cos 2 26 q = 
1 at NLO consistently which will be necessary for the cancellation of the ultraviolet 
divergences of the NLO corrections to Higgs decays into squarks. The scale of the strong 
coupling constants a s in Eqs. (fT5l [T71 [T8l [T9l [211) has been identified with the input scale 

Qo- 

In order to obtain fully consistent input parameters Eqs. ([HI EH [EH, EDD would have to 
be solved iteratively until the squark pole masses do not change anymore. Performing this 
iteration explicitly, however, shows that for certain MSSM scenarios convergence cannot 
be reached. This happens in particular for scenarios where the parameters M~ L and M| 
are close to each other so that the mixing angle 6 q is driven towards ±7r/4 and Am| 
towards —Am? . This leads to a situation in which the squark pole masses do not change 
anymore, but the mixing angles 6 q , 6 q do not correspond to the relations of Eqs. (|T3l [251) so 
that the iterated parameters are inconsistent. This situation can be avoided by stopping 
the iteration earlier before running into this inconsistent limit. Since the iteration adds 
contributions beyond NLO to the input parameters and the squark masses, the iteration 
effects can be considered as arbitrary in a NLO analysis. We have chosen a procedure 
where no iteration is performed as described above. With the shifted parameters M~ L R 
of Eq. ([20]) the squark pole masses are calculated by Eq. (fHl) with vanishing Am^ 2 
terms and the mixing angles according to Eqs. (fT3| [25]) . since the NLO corrections are 
absorbed in the shifted values for M qLR . However, for a consistent treatment of the A q 
terms of Eq. ([S]) iterations for the effective bottom mass rhb(Q ) are performed for the 
sbottom sector. This ensures that the effective bottom mass of the corrected mass matrix 
corresponds to the NLO sbottom pole masses consistently. 

The NLO neutral Higgs couplings to squarks read in the current eigenstate basis as 

dtgA^R) = rh 2 q (VR)g? + M z(hq- e q sin 2 6 w )gz 
= + M 2 z e q sin 2 6 w g% 

9tq R M = -^^[Ml-A q {»n)9t] (26) 

with the couplings gf listed in Table [2j Note that we use a common renormalization 
scale hr for the effective quark mass m q and the trilinear coupling A q for simplicity. The 
corresponding couplings to the stop/sbottom mass eigenstates gi j2 are obtained by the 
appropriate rotations according to Eq. ([7]) by the on-shell mixing angle 6 q of Eq. fl2"5l) . 
The choice of the on-shell mixing angle 6 q for these rotations implies that this is also the 
relevant mixing angle involved in the physical processes. The on-shell mixing angle is 
uniquely related to the soft SUSY-breaking input parameters according to Eqs. (I13II25P 
and thus in particular to the parameter A q . 

In our numerical analysis we include the running of the effective quark mass rh q and the 
trilinear coupling A q up to the NLL level in (SUSY-)QCD. Note that due to the counter 
terms for the effective quark masses of Eq. (1191) the heavy particles are decoupled from 
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its running. We have included 5 (6) light flavours in the NLL running of the effective 
bottom (top) masses. For the running of the trilinear couplings A q we have kept the 
contributions of all coloured particles including the top quark, the squarks and the gluino 



which is consistent with the counter term 5A q of Eq. (118]) . The running effective quark 
mass for Np = 5,6 light flavours is thus given by 

A c[a s (fi R )/TT] 



c(x) 
c(x) 



c[a s (m q )/n] 



23 
~6~ 



x 



12 

23 



' 3731 
1 + 3174 X 



7 

— x 
2 , 





137 


)' 


1 + X 




98 



for q = b 
for q — t 



5) 
6). 



(27) 



where m q denotes the corresponding quark pole mass. The NLL evolution of the trilinear 
couplings A q with Np = 6 quark and squark flavours is determined as the solution of 
the corresponding two-loop renormalization group equations (RGEs) with respect to the 
strong coupling constant [TT], [19] and can be expressed aa^| 

®s,susy(hr) 



R 



A q (Qo) + M 3 (Q ) < — — 



16 



<x s ,susy{Qq, 



- 1 



16 a S:& i si 
27 



(Qc 



7T 



a 



^ 1 a s ,susY(Qo) 

6 7T 



s^USY^R 



a 



s,SUSY 



(Qc 



with the running strong coupling a St susY including all quarks, squarks and gluinos, 



a s ,susY(fi) 



12tt 



91og(/iVA 



SUSY) 



14 log log(/i 2 / A 



susy) 



9 log(/^7AW) 



where Asusy is determined by the matching condition 



a s ,susY(Qo) = at 8 (Qo) < 1 + 



a s {Qc 



7T 



Q 



log^ + ^logS + 



6 



Ml 



-Vloj 



2 

mi 



(28) 



(29) 



(30) 



between the strong coupling a Sy susY including all quarks, squarks and gluinos and the 
low-energy QCD couplings a s with 5 active light flavours. The MS mass M 3 (Q ) of the 
gluino can be obtained from the gluino pole mass Mg, 



M 3 (Q ) = M g 



a s (M g 



3 
-( 

4 



q.t 



+ J2[Bi(Ml;m q ,m qi ) + 



r(l + e)(47r) 
2e 



1 /i 2 

+ 2 lo 4 



o 



-f-r 



Mi 



sin 29 q B (M~;m q ,m qt ) 



(31) 



8 Note that the RGEs have been transformed into the MS scheme appropriately [TT]. In the DR 
scheme the coefficient 1/6 is replaced by 7/6, the coefficient 16/27 by 4/27 and the MS gluino mass 
Ms{Qo) by the corresponding DR gluino mass in Eq. ([28]) . 
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where Ca = 3. The sum is performed over all 3 (s) quark generations. For the first 
two generations we neglected the quark masses in the squark mass matrices so that the 
corresponding squark pole masses are given by 



- 2 



Ql/r 



<„(*) + ft-"* Wn ' 



7T 



M|log^| + 



M? - m] (Qo) 



+ 2 M &/«^°) + 9 M « + ^ J IO. 



M « L/ r(Qo) 



Ml 



>(32) 



3 Results at NLO 



-< 



\ — 

s q 



Figure 2: MSSM Higgs boson decays into squark- antis quark pairs at leading order. 



At leading order the MSSM Higgs boson decays into squark-antisquark pairs are de- 
scribed by the diagram shown in Fig. [2J The corresponding partial decay width reads 
as 

N C G F 



with the two-body phase-space function 



(33) 



Ay 



2 

Ml 



rn 



Ml 



Ml 



(34) 



While Higgs boson decays into squarks of the first two generations are generally sup- 
pressed, Higgs decays into sbottoms and stops develop sizable branching ratios, whenever 
they are kinematically possible [7] . In particular scenarios their branching ratios can reach 
a level of 80-90% [20]. 

The SUSY-QCD corrections at NLO have been calculated in Refs. [8]. However, the 
analytical results will be repeated in this section in our notation and for the investiga- 
tion of calculational details once consistent input parameters at NLO are included in the 
theoretical analysis. The SUSY-QCD corrections consist of one-loop virtual corrections, 
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Figure 3: Virtual corrections to MSSM Higgs boson decays into s quark- antis quark pairs 
at next-to-leading order. The external self-energy blobs represent the self-energy diagrams 
ofFzg.m 



mediated by the Feynman diagrams depicted in Fig. [3j supplemented by the renormal- 
ization of the mass, coupling and wave function parameters involved, and real corrections 
due to gluon bremsstrahlung, see Fig. @J We have performed the loop and phase-space 
integration within dimensional regularization in n = 4 — 2e dimensions so that ultra- 
violet and infrared singularities appear as poles in e. The infrared singularities cancel 
after adding the real and virtual corrections. The ultraviolet poles disappear after adding 
the corresponding counter terms which will be discussed in detail in the following, since 
they have to correspond to the schemes in which the squark masses, mixing angles and 
couplings are defined. 

The wave-function counter terms 5Za (i = 1,2) are determined by normalizing the 
residues of the diagonalized squark propagators to unity. The explicit calculation of the 
diagrams of Fig. [T] leads to the expressions 

6Z U = {S (m|; 0, nig.) - S (m|; M~ g , m q ) + 2m\ff Q {m\\ 0, m & ) 

+ [M 2 g +m 2 q -ml + (-iy2M g m q sm26 q ] B' {m 2 q - M~ g ,m q )} , (35) 

where the derivative of the two-point function is defined as 

d 

B' (p 2 ; m 1 ,m 2 ) = —B (p 2 ; m 1 ,m 2 ) ■ (36) 
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Figure 4: i?eaZ corrections to MSSM Higgs boson decays into s quark- antis quark pairs at 
next-to-leading order. 



Throughout our calculation in this section the quark mass m q is defined as ^(/i^) of 
Eq. OH]) for decays into sbottom pairs and as m t (/i#) of Eq. fl9]) for decays into stops. 
The mixing angle is renormalized by the anti-Hermitian counter term of Eq. ( I23p . The 



trilinear couplings A q (fi R ) (q = b,t) are defined in the MS scheme at the scale which 
is identical to the running couplings in the DR definition at NLO. Thus, the counter term 
of A q (fj,ji) can be cast into the form 



5A q = C F -r(l + e)(4vr) £ M § {- +log %■ \ . 

71 I e ) 



(37) 



The quark mass counter term 5m„ is given by 



5m„ = Srh n 



(38) 



with the top and bottom mass counter terms 5rh q of Eq. (fl9l) with Q = fiR. 

The calculation of the diagrams of Fig. [3] for the virtual corrections and those of Fig. 0] 
for the real corrections leads to the final result after renormalization, 
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B (m qi ; 0, m 9 -J + B (rrig. ; 0, m^) - B (M^; m qv m 
-2(M| - m| - m|.)C (m|, M|, mf ; 0, m ffi , m^) 



B (rn? ; M §) m ? ) + B (mhM~ g , m q ) + 25 (Mj^; m„ m g ) 



+(2M| + 2mJ - m| - m|)C (m|, A^, m| .; M 9 , m q , m q ) 
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m q sin 20 q 
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\M 2 
1 



^ff^. - 

A-m 4- A, 4 

0ro„ ^ + 0A 6Aq 



■q ^^q 
1 - Pi - Pi 



2/3 



logXo PiPjXo 

h log x log ■ 



--log xx - - log x 2 + ALi 2 



-x 



, 1 , 1 APi , , 1 + A + Pi , 



/? 4 

- 2Li 2 (l - xi) - 2Li 2 (l - x 2 ) 

pifogXi + pj log x 2 



2/3 



/3 



(39) 



where in C rea i we used the abbreviations /3 = ^^/Xij with the two-body phase-space func 
tion Ajj of Eq. ( 1511) and 



;r 



1 - Pi ~ Pj - P 



Xi 



1 - pi + pj- /3 



x 2 



1 + Pi ~ Pj ~ P 



(40) 



1 — Pi — Pj+P 1 - pi + pj + (3 " 1 + pi - pj + f3 

with pi = m 2 ~/M^ > . The two mixing matrices 1Z and S used in the coefficients C 2 ' H and 
C 3 are given by 



ft 



sin 26L 



cos 26L 



cos 2# g — sin 26 q 



S 



COS 29 q 

— sin 2#„ 



• sin 20 q 
cos 26*0 



The three-point function Cq is defined as [16 

Co(pf,P 2 i ,Pn;rn 1 ,m 2 ,rn3) 



-i(47r) 2 /i 2£ 



(2vr) n [fc 2 - m 2 ] [(jfc + Pl ) 2 - m 2 ] [(jfc + p 12 ) 2 - m\ 



(41) 



(42) 



with P12 = pi + p 2 . The finite remainders of the mixing angle renormalization at the 
external squark legs can be cast into the form of the matrix 



T 



9(7* 4- (7* 

y <Jl<?2 y <?l<?l 1 y <?2 

?_ 4-0?- 2a?- 



(43) 
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multiplied by the finite shift 



A9 q = , f - 69 q = -SRe 121 " 2 , , 44 



which includes the corresponding ant i- Hermit ian counter term 59 q of the mixing angle as 
given in Eq. f[2"3l . In Eq. (|44l) the singularity for m| 2 — )■ m? x is cancelled by the anti- 
Hermitian counter term. Finally the derivatives of the squark couplings gf iq -. to the Higgs 
bosons are given by 

do h '? h,H 
-^L = 2m q g h q H 5 l3 + 9 -^n i3 
om q m q 



da-: 11 



dA q 2 » q 13 

tyqiqz ^9q 2 q\ _ 9q\q 2 



dm q dm q m q 

d 9iq 2 d 9i qi _m q . 



^9? ■ (45) 



dA q dA q 2 

The results for the coefficients C* are ultraviolet and infrared finite after all individ- 
ual contributions are added up since our renormalized mixing angles fulfill the relation 
sin 2 29 q + cos 2 29 q = 1 consistently. The scale of a s in Eq. (I39j) is identified with the renor- 
malization scale everywhere apart from the term within the effective bottom mass of 
Eq. OH]). It should be noted that our final results do not contain any leading tg/3-enhanced 
corrections of the A^-type anymore due to our scheme choices of the effective masses and 
couplings at NLO. Due to this property the problems of a naive renormalization program 
have been solved consistently leading to moderate radiative corrections to the Higgs decay 
widths into squark-antisquark pairs as will be shown explicitly in the next section. 

By using power counting arguments along the lines of Ref. [13] it can be shown that 
our approach resums all leading terms of 0(a"/z n tg n /3) so that the residual HO corrections 
are free of these leading terms. In the current eigenstate basis there are three sources of 
terms proportional to tg/3: (i) the Higgs couplings to squarks and quarks as in Table Q] 
and Eq. ( 126]) which drop out in the relative corrections; (ii) off-diagonal mass insertions in 
the sbottom propagators which are power-suppressed by mb/Ms due to the KLN theorem 
and the m^/Ms structure of these insertions; (iii) potential tg/3-enhancements of the 
counter terms. The MS renormalization of as in Eq. ( 1371) does not involve any tg/3 
enhanced corrections as can also be inferred from the beta functions of the RGEs for the 
trilinear couplings up to three- loop order [T9] . The tg/3 enhancement of the bottom-mass 
counter term entering the corresponding bottom Yukawa couplings involved in the Higgs 
couplings to bottom quarks and the sbottom mass matrix is explicitly absorbed in the 
effective bottom mass of Eq. (jSJ) which contains the corresponding threshold correction 
A;, in resummed form [12] [13] . Therefore the counter terms Srrib and 5Ab are free of tg/3 
enhanced corrections at all orders up to terms which are suppressed by nib/ Ms- 
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Close to the threshold M$ 
Coulomb singularities, 



+ rriq. the NLO SUSY-QCD corrections develop 



7T~ 



Ay {rnq.+mq^ 



for M$ — >■ m„- + nij 



«3 ' 



(46) 



which agrees with the usual Sommerfeld rescattering correction factor [21] . These singu- 
larities for a/A^ — y can be regularized by taking into account the finite decay widths 
of the squarks. Moreover, these Coulomb factors can be resummed systematically. Both 
effects are expected to be relevant close to threshold, but are beyond the scope of this 
paper. 

On the other hand far above the threshold M$ 3> m^. ^> m q the NLO corrections 
approach the asymptotic limit 



dg, 



4> 



+2 



m s daU 

9k dA « 



9 9t 



- 1 



m q cos 26 q 
+ U] M, 



<S 9, 







3, ul 1 Ml 7 
4 l0g M|~4 l0g M| + I + 79 



L s 



g M 2 S + 4 g M 2 S + 2 



log 



Ml 



- 2 



+ 



h,H "V 

9Mj 



Sij log 



Ml 



, 1 , M i 

1 - 7T log 



+ 5d 



(47) 



where we have identified all supersymmetric masses, Ms = m, 
glected the quark mass against the other masses. The term 7 g is given by 



Mg, and ne- 



7t 



A - /i/tg/3 
4Mc 



A, 



7fe 



4M S 



(4E 



thus confirming the absence of large corrections of 0(a s /itg/3jj- The asymptotic result 
far above the threshold develops a double logarithmic contribution in the second square 
bracket of Eq. (T47j) which originates from the second diagram of Fig. [3] in the Sudakov 
limit. It cannot be absorbed in an appropriate choice of the renormalization scale fiR but 
will cancel against the corresponding double logarithm of gluino radiation in association 
with quark-squark pairs in the large Higgs mass limit which has not been taken into 
account in this work. In the following we will identify fi^ with the corresponding Higgs 
mass Mq, as the central scale choice. 



9 In the on-shell scheme, which relates the renormalization of A q to the on-shell quark and squark mass 
as well as the mixing angle counter terms as in Eq. (|21l) and renormalizes the quark mass on-shell, large 
corrections of C(a s /z 2 tg 2 /3) emerge for the partial Higgs decay widths into squarks. 
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4 Numerical Results 



The numerical analysis of the neutral Higgs boson decays into stop and sbottom pairs is 
performed for two MSSM scenarios, one close to the one of Ref. [10] where we lifted the 
gluino and squark masses of the first two generations beyond the mass bounds from the 
LHC [22], and a second scenario with large SUSY-breaking masses and large mixing in 
the stop sector as representative cases: 

A) Qo = 300 GeV, 

tg/3 = 30, M h (Q ) = M~ bL (Qo) = M %R (Q ) = 300 GeV, 

M iR (Q ) = 270 GeV, M bR (Q ) = 330 GeV, 

M^ R (Qo) = M 1lr (Q ) = M~ g = 1 TeV, A t (Q ) = 150 GeV, 

A b (Q ) = -700 GeV, A T (Q ) = 1 TeV, M 2 = 1 TeV, y, = 260 GeV 

B) Qo = 500 GeV, 

tg/3 = 30, M h JQo) = M- bLR (Q ) = MfUQo) = 500 GeV, 
Mq L>R (Qo) = M iL R (Q ) = Mg = l TeV, A t (Q ) = A b (Q ) = -1.5 TeV, 
A T (Q ) = 0, M 2 = 500 GeV, y = 500 GeV , (49) 

where Mq L R (Q ) and Mi lr {Qq) denote the squark and slepton mass parameters of the 
first two generations. The results of this work have been implemented in the program 
HDECAY [23] which calculates the MSSM Higgs decay widths and branching ratios in- 
cluding the relevant higher-order corrections |20j. We use the RG- improved two- loop 
expressions for the Higgs masses and couplings of Ref. [23] which yield predictions for 
the Higgs boson masses that agree with the diagrammatic calculations of Ref. [2] within 
3-4% in general. Thus the leading one- and two-loop corrections have been included in 
the Higgs masses and the effective mixing angle a. Consistency of our scheme with the 
scheme and scale choices of Ref. [23] requires the evolution of our sbottom parameters A b 
and Mr to the scale 

°L/R 



Q b = max 




and the stop parameters A t and to the scale 

Qt = max ( ^M 2 h {Q Q ) + m 2 t {m t ), \Jm] r {Q q ) + m 2 t {m t )\ . (51) 
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Since the scales Q and Q b / t are of similar order of magnitude we neglect resummation 
effects so that the relations are given b}q 



A b {Q b ) 
M~ bR {Q h ) 

MQt) 

M 2 iL (Q t ) 
Ml(Q t ) 
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with the abbreviations 

a b 

a t 
X b 



m 2 b (Qc 



2ttv 2 cos 2 (3 

rn 2 t(Qo) 
2ttv 2 sin 2 [3 



M h (Q ) + M~ bR (Q Q ) + M 2 Hi + Al(Qo) 
M 2 iL (Q ) + Ml(Qo) + Ml 2 + A 2 t (Q ) 

and the soft SUSY-breaking Higgs mass parameters 



M 2 Hl 



M|sin 2 /3 



Ml 



cos 2(3 — // 2 



M 



H 2 



M\ cos 2 /3 + cos 2 /3 



A* 



Q* 



Qo 



loS Ql 



(52) 



(53) 



(54) 



The top pole mass has been taken as m t = 172.6 GeV, while the bottom quark pole 
mass has been chosen to be m b = 4.60 GeV, which corresponds to a MS mass m b {Tn b ) = 
4.26 GeV. The strong coupling constant has been normalized to a s (Mz) — 0.118 which 
corresponds to the QCD scale A 5 = 226.2 MeV for 5 active flavours. The related SUSY- 
QCD scale Asusy of Eq. ( |29l) amounts to Asusy = 2.028 keV in scenario A and Asusy = 
1.756 keV in scenario B. The MS masses involved in the scales Qt, b can be derived as 
m t (m t ) = 165.1 GeV and m b {m t ) = 2.78 GeV while the MS gluino mass at the input 
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Note that the left-handed soft supersymmetry-breaking squark mass parameters M b (Qb) and 
Mi L (Qt) are not equal, because they are evaluated for different scales Q t and Qb- We have modified the 
calculation of Ref. 24 to account for these differences consistently. 
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scale Q amounts to M 3 (Q ) = 1.028 TeV in scenario A and M 3 (Q ) = 1.006 TeV in 
scenario B. The squark masses in the two scenarios amount in particular to 



A) m h = 165.5 GeV, m h = 278.0 GeV, m~ bi = 182.5 GeV, m~ h% = 265.6 GeV, 
mu L = 989.8 GeV, m^ R = 990.7 GeV, m s = 993.1 GeV, m~ dR = 991.6 GeV 



B) m h = 202.1 GeV, m h = 708.0 GeV, m~ hi = 473.0 GeV, m h = 536.3 GeV, 

rriu L = 1.011 TeV, m^ R = 1.012 TeV, = 1.014 TeV, m~ dR = 1.012 TeV , (55) 

where fn i i L/R and m d L/R denote the up- and down-type squark masses of the first two 
generations. 

In Fig. [5] we display the results for the partial decay widths of the heavy neutral 
Higgs particles into stop pairs and the relative corrections in scenario A. Since the stops 
are moderately heavy all decay modes are kinematically allowed for Higgs masses above 
about 560 GeV. The partial decay widths range at the few GeV level and turn out to be 
similar for all Higgs bosons above the corresponding thresholds except the scalar Higgs 
decay widths into non-diagonal stop pairs which are smaller. Apart from the Coulomb 
singularities at threshold the SUSY-QCD corrections are of moderate size in our scheme. 
Similar results have been obtained for the neutral Higgs boson decay widths into sbottom 
pairs as shown in Fig. [6J The moderate size of the full SUSY-QCD corrections confirms 
the proper absorption of the A& terms as in Eq. (jSj). If the trilinear coupling Af, and 
the bottom quark mass would be renormalized in the on-shell scheme the SUSY- 
QCD corrections would increase the LO results by more than an order of magnitude so 
that the result in the on-shell scheme becomes totally unreliable [TUl ES]- The partial 

widths of the non-diagonal Higgs decays H — > bib 2 , b 2 bi are small, since the mixing factor 
cos2# fe is small and thus the contributions of the left-right couplings g?- are suppressed, 

while the diagonal couplings q^ 1 z and q? r cancel each other to a large extent in the 

o x- o ^b L b L ^b R b R 

evaluation of the coupling g~*~ of Eq. (J7J). Using the scheme of Ref. [17] for the mixing 
angle, i.e. defining the mixing angle counter terms as 56 q = ^eE^fQ 2 ) /jVn| 2 ~ m |j) with 
Q 2 = m| x or Q 2 = m~ 2 , the results agree with ours within less than l%rH 

In order to obtain an estimate of the residual theoretical uncertainties we show the 
renormalization scale dependence of the partial decay widths into stop and sbottom pairs 
in Fig. [7J The scale dependences of the partial decay widths into stop and sbottom pairs 
are significantly reduced at NLO. Varying the renormalization scale by a factor of 2 around 
the central scale [Mr = M h /a the theoretical uncertainties can be estimated as 5-10% for 
decays into stop and sbottom pairs at NLO. 

Fig. [8] displays the corresponding branching ratios of the dominant decay modes of 
the heavy scalar and the pseudoscalar Higgs boson in scenario A. It is clearly visible that 
the decays into stop and sbottom pairs belong to the dominant decay modes for masses 
above about 300-400 GeV, i.e. where they are kinematically allowed, reaching branching 



The results using the tree- level like mixing angle q of Eq. ([13]) develop sizeable differences to the 
results in our scheme for squark masses 2 close to each other, where the partial decay widths H — > 
<?i<?2, <72<7i turn out to be negative thus signalizing a basic problem with the tree-level-like mixing angle. 
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Figure 5: SUSY-QCD corrected partial decay widths (upper) and the relative SUSY-QCD 
corrections (lower) of the heavy scalar and the pseudoscalar MSSM Higgs boson decays to 
stop pairs as functions of the corresponding Higgs masses for tg/3 = 30 in scenario A. The 
full curves include the SUSY-QCD corrections, while the dashed ones are the LO results. 
The small kinks originate from the thresholds of the other stops in the virtual corrections. 
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Figure 6: SUSY-QCD corrected partial decay widths (upper) and the relative SUSY-QCD 
corrections (lower) of the heavy scalar and the pseudoscalar MSSM Higgs boson decays 
to sbottom pairs as functions of the corresponding Higgs masses for tg/3 = 30 in scenario 
A. The full curves include the SUSY-QCD corrections, while the dashed ones are the LO 
results. The small kinks originate from the thresholds of the other sbottoms in the virtual 
corrections. 
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Figure 7: Scale dependences of the partial decay widths into stop pairs (upper) and sbottom 
pairs (lower) as functions of the renormalization scale in units of the central scale /i = 
Mh/a i n scenario A. The broken lines show the leading- order scale dependences and the 
full curves the NLO results. 
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ratios of up to about 50% in total. Moreover, in scenario A the decay modes into f's and 
neutralinos play a significant role, too. 

In Fig. IH] we display the final results for the scenario B. Since the masses m^ 2 ,m^ i2 

are close to or larger than 500 GeV, the decay modes H — y t 2 t 2 , b\b 2l b 2 b\, b 2 b 2 are kine- 
matically forbidden for Mn < 1 TeV. The partial decay widths of decays into stop and 
sbottom pairs range at the few GeV level as in scenario A once they are kinematically 
allowed. It can clearly be inferred from Fig. [9] that the SUSY-QCD corrections are of 
moderate size apart from the threshold regions where the Coulomb singularities enhance 
the size of the corrections. It should be noted that the results with the sbottom mixing 
angle renormalized via the tree-level relation of Eq. (^T§ would lead to unphysical nega- 
tive partial decay widths for heavy scalar Higgs decays H — y bib 2) b 2 b\ for Higgs masses 
larger than 1 TeV, where these decay channels open up, so that this scheme is strongly 
disfavoured. These negative contributions can be traced back to the uncancelled artificial 
singularity for rag ~ mg in the finite remainders A6b of Eq. (|4*4"]) . if 58b is replaced by 
the counter term 56 b of Eq. f l2T]) [ID]. 

In Fig. [TOj we show the branching ratios of the heavy scalar and the pseudoscalar 
Higgs boson in scenario B as functions of the corresponding Higgs masses. As in scenario 
A the decay modes into stop and sbottom pairs play a significant role, once they are 
kinematically allowed. This starts to be the case for the heavy scalar Higgs boson H 
already for masses above about 400 GeV, where the decay into light stop mass eigenstates 
opens up. The pseudoscalar Higgs boson can only decay into mixed pairs of a light and 
heavy stop and sbottom eigenstates so that these decay modes open up for masses above 
900 GeV. The decays into charginos and neutralinos contribute significantly for larger 
Higgs masses. 

For smaller values of tg/3 the decay modes H — y hh, ti and A — y Zh, ti play a significant 
role, while the decays into sbottom pairs are usually suppressed. However, the neutral 
Higgs boson decays into stop pairs still play a significant role and can even be the dominant 
heavy Higgs boson decays in a large Higgs mass range. The SUSY-QCD corrections to 
these decay modes are of similar size as in the scenarios with large values of tg/3 analyzed 
in this work. 

5 Conclusions 

In this work we have analyzed the NLO SUSY-QCD corrections to heavy MSSM Higgs 
decays into stop and sbottom pairs and elaborated a consistent scheme and setup for the 
theoretical NLO results and the corresponding input parameters, i.e. the squark masses 
and their couplings to the Higgs bosons. A reliable determination of these corrections 
requires a proper treatment of contributions which are enhanced by tg/3. The derivation 
of the input squark masses and couplings with NLO accuracy requires a modification of 
the tree-level relations among the stop and sbottom masses by higher-order corrections, if 
they should correspond to the physical pole masses. Since the size of these modifications 
depends on the MSSM parameters and pole masses themselves an iterative procedure 
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Figure 8: Branching ratios of the heavy scalar H and the pseudoscalar A MSSM Higgs 
bosons in scenario A as functions of the corresponding Higgs masses. The curves for 
decays into neutralinos x°X° an d squarks qq represent the corresponding sums over all 
possible mass eigenstates. 
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Figure 9: SUSY-QCD corrected partial decay widths (upper) and the relative SUSY-QCD 
corrections (lower) of the heavy scalar and the pseudoscalar MSSM Higgs boson decays 
to stop and sbottom pairs as functions of the corresponding Higgs masses for tg/3 = 30 
in scenario B. The full curves include the SUSY-QCD corrections, while the dashed ones 
are the LO results. The small kinks originate from the thresholds of the other stops in the 
virtual corrections. 
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Figure 10: Branching ratios of the heavy scalar H and the pseudoscalar A MSSM Higgs 
bosons in scenario B as functions of the corresponding Higgs masses. The curves for 
decays into charginos x + X~ > neutralinos x°X° an< ^ squarks qq represent the corresponding 
sums over all possible mass eigenstates. 
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should be applied. These iterations, however, do not converge for all scenarios so that 
we restricted the evaluation to a single iteration step. The final setup of the input pa- 
rameters matches NLO accuracy consistently. A preferred choice of the scheme for the 
stop and sbottom mixing angles is the on-shell definition via the ant i- Hermit ian counter 
term which removes spurious singularities for nearly mass degenerate cases which drive 
the off-diagonal decays — ^ 6162 to negative and thus unphysical values. The final NLO 
corrections are of moderate size apart from the expected Coulomb singularities at the 
decay thresholds. The residual theoretical uncertainties at NLO have been estimated to 
be 5-10%. 

Note added in proof. During completion of this work Ref. [26] appeared which dis- 
cusses several options for the renormalization of the stop and sbottom sector with NLO 
accuracy. Their treatment of the stop and sbottom sectors is different from our approach 
in particular for the effective bottom mass and the squark mixing angles. 
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